Reaction-diffusion processes on correlated and uncorrelated scale-free networks 
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We compare reaction-diffusion processes of the A + A^O type on scale- free networks created with 
either the configuration model or the uncorrelated configuration model. We show via simulations 
that except for the difference in the behavior of the two models, different results are observed within 
the same model when the minimum number of connections for a node varies from fc m i n = 1 to 
fcmin = 2. This difference is attributed to the varying local properties of the two systems. In all 
cases we are able to identify a power law behavior of the density decay with time with an exponent 
/ > 1, considerably larger than its lattice counterpart. 

PACS numbers: 82.20.-w, 05.40.-a, 89. 75. Da, 89. 75. He 
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Recently, a growing interest for dynamic processes tak- 
ing place on scale- free networks has arisen jj], 0] ■ A scale- 
free network is one where the node connectivity k (num- 
ber of links on a node) follows a power-law distribution 
of the form 



P(k) 



(1) 



where 7 is a positive number, typically in the range 
2 < 7 < 4. In this frame, we recently presented Q 
simulation results for reaction-diffusion processes both 
of the type A + A — > and A + B — ■> 0, where the sub- 
strate for diffusion is a scale-free network. Following this, 
Catanzaro et al. Q developed an elegant theory for the 
A + A — > process, which applies on networks created 
with the uncorrelated configuration model (UCM) 
Their analytic results were found to be in good agree- 
ment to their simulations, but were deviating from the 
results in Ref . Q , where the networks were created with 
the configuration model (CM). The authors attributed 
the difference of the results on the different method of 
preparing the networks. In this Brief Report we directly 
compare results for both the UCM and CM models. 

The Configuration Model (CM) has become more or 
less the standard for simulating networks of a given 7 
value in the literature. For each one of the N system 
nodes a random k value is assigned, drawn from the prob- 
ability distribution function of Eq. Q . Pairs of links are 
then randomly chosen between the nodes, taking care 
that no double links between two nodes or self-links in the 
same node are established. When the maximum value of 
k is not pre-defined, then the natural upper cutoff scales 
with the system size as fc max ~ jyVa- 1 ). This model 
is known to create correlations in the connectivity dis- 
tribution between the system nodes for 7 < 3, in the 
sense that the average connectivity of a node's neighbors 
depends on its rank k. In other words, nodes that are 
highly connected prefer to attach to nodes with lower k- 
values, rather than to equally well-connected nodes. The 
Uncorrelated Configuration Model (UCM) uses the same 
construction algorithm, with the difference that the up- 
per cutoff is fixed in advance to fe max = N 1 / 2 . Then, it 
has been shown |f| that connectivity correlations dissa- 
pear and the average connectivity for the neighbors of 



any node is constant. 

In Ref. we had found that the reaction rate was 
evolving surprisingly faster than on regular lattices. The 
variation of the particle concentration p(t) was still fol- 
lowing a power law with time t of the form 
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but with a value / > 1 for 7 < 3.5. This result was 
valid asymptotically within the simulation accuracy, and 
before finite-size effects settled in. The networks that 
we had used had been created with the CM model with 
a minimum value for the connectivity of a node (lower 
cutoff) fc ro i n = 1. Since this process may create isolated 
clusters, we only had used the largest cluster which, de- 
pending on the value of 7, would span from 35% to 100% 
of the system nodes. 

In Ref. H| the A + A — > process was studied, and it 
was analytically found that for an infinite size network 
the exponent / of Eq. J2J is given by / = 1/(7 — 2) when 
2 < 7 < 3, i.e. again / > 1. For 7 = 3.0 the behavior was 
predicted analytically to be 1/p ~ tint. However, for fi- 
nite size networks and long times the behavior of 1/p is 
masked for all 7 by the mean- field exponent / = 1, which 
seems to also be valid asymptotically from the simulation 
results. Both the analytic solution and the simulations 
in that paper were based on networks created under the 
UCM model, with a lower cutoff value of fc m i n = 2. The 
authors of Ref. J4J reported that in their simulations it 
was not possible to find a noticeable regime with / > 1, 
and they attribute the discrepancy between the two stud- 
ies solely in the different network creation method. We 
here argue that a power law regime with / > 1 is indeed 
present and clearly identified in all cases. Additionally, 
a very important factor is the lower cutoff value, fc m ; n - 

In the present Report we present and compare simu- 
lation results for both the CM and UCM models with 
fcmin = 1 and fcmin = 2. Results for the time evolu- 
tion of the particle density in all four possible cases are 
shown in Fig. ^ All four curves in the case of 7 = 2.5 
(Fig. [TJi) behave differently from each other. The CM 
and UCM models clearly yield different results, but even 
within UCM or CM the curves for fc m ;. 
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FIG. 1: Plot of the reaction progress 1/p — 1/po, as a function 
of time for the A+A— >0 reaction on scale-free networks of (a) 
7 = 2.5, and (b) 7 = 3.0. Results correspond to CM (thin 
lines) and UCM (thick lines) models, with fc m i n = 1 (solid 
lines) and fc m in = 2 (dashed lines). The initial density was 
po = 0.5. All results correspond to networks of TV = 10 6 
nodes. 
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FIG. 2: Results for networks with 7 = 2.5. Plots of inverse 
particles density (1/p — 1/po) divided by (a) t and (b) t', as 
a function of time. Thick lines correspond to the UCM and 
thin lines to the CM model while in solid curves k m - m = 1 and 
in dashed lines fc m in = 2. In (b) the / values are shown on 
the plot. 



deviate from each other. The main difference is that the 
crossover to the power law behavior for fc m ; n = 1 appears 
later in time, and the asymptotic mean-field behavior 
(/ = 1) also exhibits itself roughly one decade later. 

The case of 7 = 3.0 (Fig. is simpler. For this 
7 value the two models (CM and UCM) are expected 
to coincide, since in general the natural upper cutoff in 
the CM model scales as k c ~ TV 1 ' 7-1 , and for 7 = 3.0 
the value k c = N 1 / 2 is exactly the same as in the UCM 
model. This coincidence is shown to be valid in the figure. 
However, the results for fc min = 1 are still different than 
the results for fc m ; n = 2. The crossover to the power-law 
behavior is not as prominent as in the case of fc m j n = 1. 

We then assess whether there exists an observable 
power-law behavior in these curves and whether asymp- 
totically a linear regime masks this power-law. The re- 
sults for 7 = 2.5 are presented in Fig. [21 In the first 
plot we divide the curves of Fig. ^ by the linear func- 



tion t. If asymptotically the expected behavior is linear 
(/ = 1 in Eq. [2J then we expect that this division will 
yield a constant value, i.e. a horizontal line parallel to 
the x-axis. The only case which is close to this behavior 
is that of the UCM model and fc m i n = 2 (thick dashed 
line in the plot). Even then, the result is not entirely 
satisfactory, and the curve seems to have a slope greater 
than 0. Notice also that although from Fig. one can 
visually get the impression that a slope of 1 is describing 
quite accurately this curve, the more detailed analysis in 
Fig. [2^ reveals that this is not true. In all other cases, 
there are only weak hints for a parallel line, but it cannot 
be argued that this linear relation is in general valid, and 
this is why it was not observed in Ref. [3J. The abruptly 
falling curves at longer times are due to very low parti- 
cle densities, where less than 10 particles remain in the 
system. 

In Fig. |2t>, we also test the hypothesis for power-law 
behavior by dividing the raw data of Fig.QJiby t f . bince 
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FIG. 3: Results for networks with 7 = 3.0. Plots of inverse 
particle density (1/p — 1/po) divided by (a) and (b) tint, 
as a function of time. The CM and UCM models yield the 
same curves, so we only present results for CM with fc m i n = 1 
(solid lines) and fc m in = 2 (dashed lines). 

curves corresponding to different models exhibit varying 
exponents / we choose the value of / that maximizes the 
extent of the horizontal part in each line. 

The presence of a power law regime is clear in all oc- 
casions, but it is also noteworthy that the validity of the 
power law lies in a narrower time range for A: m i n = 2, and 
especially from the UCM model, rendering this estima- 
tion more difficult. Most probably, this is the reason that 
the authors of Ref. ^4J did not characterize this regime. 
The case where in a simulation a particular power law 
shows only in a very limited time domain is not unusual, 



see e.g. the Zeldovich behavior of the A + B reaction on a 
3-dimensional lattice, where the density exponent has a 
nominal value of / = 0.75, which shows only in one time 
decade |6J . The slopes of the curves for A: m i n = 1 are con- 
sistently larger compared to the slopes of the fc m i n = 2 
curves. We also note that the analytical prediction of 
Ref. Q predicts a slope of / = 2.0, which is very close 
to the slope of the CM model for A: m i n = 1, but deviates 
remarkably from the slope of the UCM model, although 
this prediction was based on the uncorrelated network 
hypothesis. 

We have already seen that when 7 = 3.0 the CM 
and UCM models are identical. For clarity in Fig. [3] we 
analyze only the results of CM. Again, the hypothesis 
that asymptotically there is a linear increase is not di- 
rectly supported from our simulation results. On the 
contrary, if we consider a power law fit for the intermedi- 
ate time regime, we observe an extended horizontal part 
that spans a few decades. The slope for both fc m i n = 1 
and 2 are similar and close to 1.3. An increase in the 
network size N has been shown to further extend the re- 
gion of validity for the power-law behavior in the / > 1 
regime (see e.g. Fig. 1 of and Fig. 6 of 3). 

The analytic considerations in Ref. Q for 7 = 3.0 pre- 
dict a logarithmic correction (1/p ~ tint). In general, it 
is not easy to distinguish between this form and a power 
law curve. However, when we divide our raw simula- 
tion data with the tint function in Fig. Et>, we observe 
that for fc m in = 2 there is an excellent agreement, but this 
function fails to reliably describe the results for fc m i„ = 1 . 

The most probable explanation for the difference in the 
results for the varying fc m j n values is that the local envi- 
ronment in the case of fc m ; n = 1 is remarkably different 
than when k m i n — 2. Although fc m j n = 2 is known to 
yield one giant connected cluster [jj, its structural char- 
acteristics seem to be modified, basically because of the 
tree-like features (many nodes are connected to just one 
neighbor). Bottlenecks and revisitations are, thus, more 
frequent than in the case where a larger lower cutoff is 
used. 

In short, the value of the exponent characterizing a 
reaction-diffusion process on scale-free networks depends 
not only on the network construction model, but is also 
sensitive to the lower cutoff value for the connectivity 
&min, although such a dependence has not been predicted 
analytically. More importantly, we have verified the exis- 
tence of a power law regime that is large enough to be ob- 
servable in all cases, while a careful analysis also revealed 
that an asymptotic linear behavior is not in general valid. 
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